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High precision vector magnetometry with uniaxial quantum centers in silicon carbide 
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We show that uniaxial color centers in silicon carbide with hexagonal lattice structure can be 
used to measure not only the strength but also the polar angle of the external magnetic field with 
respect to the defect axis with high precision. The method is based on the optical detection of 
multiple spin resonances in the silicon vacancy defect with quadruplet ground state. We achieve 
a perfect agreement between the experimental and calculated spin resonance spectra without any 
fitting parameters, providing angle resolution of a few degrees in the magnetic field range up to 
several millitesla. Our approach is suitable for ensembles as well as for single spin-3/2 color centers, 
allowing for vector magnetometry on the nanoscale at ambient conditions. 


PACS numbers: 76.30.Mi, 71.70.Ej, 76.70.Hb, 61.72.jd 

Optically addressable atomic-scale spin centers consti¬ 
tute the basis for nanomagnetometry with high sensitiv¬ 
ity and high spatial resolution mm- The most promi¬ 
nent example is the nitrogen-vacancy (NV) defect in dia¬ 
mond and several benchmark experiments have been per¬ 
formed using this system ms, including proton nuclear 
magnetic resonance on the nanometer scale [HIT]. The 
principle of magnetometry with spin-carrying color cen¬ 
ters is based on optical detection of magnetic resonance 
(ODMR), subject to external magnetic field. In case of 
individual NV defects with spin S' = 1, the projection of 
the magnetic field on the defect axis is measured. The NV 
defect in the diamond cubic lattice is oriented along one 
out of four (111) crystallographic axes and, therefore, us¬ 
ing ensemble experiments the magnetic field vector B can 
be reconstructed mis]. Ensembles of the NV defects are 
also suggested for the implementation of high precision 
magnetic field sensors with femtotesla sensitivity [iQiin] 
and solid-state frequency standards [12]. These imple¬ 
mentations require high homogeneity of the NV centers. 
The NV defects can be fabricated with preferential align¬ 
ment [T51 E], and using nonlinear shift of the ODMR 
lines in relatively high magnetic fields of several tens of 
millitesla the transverse field component can be recon¬ 
structed m- However, in many demanding applications 
much lower magnetic fields should be detected, and the 
information on the magnetic field orientation is difficult 
to extract in this approach. 

Here, we demonstrate an alternative approach to im¬ 
plement vector magnetometry for magnetic fields below 
several millitesla, which is suitable for ensemble as well 
as for individual uniaxial spin centers with S = 3/2 [T6|. 

As a model system, we consider a silicon vacancy (Vsi) in 
silicon carbide (SiG) pTHSO] . Due to the polymorphism 
of SiG, there is a large variety of vacancy-related defects 
with appealing quantum properties nniiiDEi- All ex¬ 
periments presented here have been performed at room 
temperature on a 4H-SiG bulk crystal, possessing hexag¬ 


onal lattice structure. The crystal has been grown by 
the standard sublimation technique, such that the [0001] 
crystallographic direction (c-axis) is inclined at an angle 
of 7° to the surface normal, i.e., to the z-axis of the lab¬ 
oratory coordinate system [Fig. [^a)]. The crystal has 
been irradiated with neutrons (5MeV) with a fluence of 
10^^cm“^, in order to generate Vsi defects. Their pres¬ 
ence is verified by the characteristic photoluminescence 
(PL) in the near infrared spectral range [34] . 

The Vsi defects in hexagonal SiG have spin-3/2 ground 
state, which is split in two spin sublevels ms = ±1/2 and 
ms = ±3/2 at zero magnetic field [27]. The zero-field 
splitting 2D between these sublevels depends on the lat¬ 
tice site and polytype. In 4H-SiG there are two nonequiv¬ 
alent lattice cites and hence two different Vsi defects. 
They are distinguished by their spectroscopic features 
and labeled as VI and V2 centers [18]. Here, we present 
results for the Vsi(V2) center with 2D/h = 70 MHz 
[Fig. [^a)]. A detailed characterization of the system is 
presented elsewhere [UlSSlEilEIl [32 [34] 

A diode laser operating at 785 nm is used to opti¬ 
cally pump the Vsi(V2) center in the host SiG with 
a bandgap of 3.23 eV. The optical excitation followed 
by spin-dependent recombination results in preferential 
population of the ms = ±1/2 sublevels [35] • The PL 
from these centers is passed through 800-nm and 850-nm 
longpass filters and detected by a Si photodiode (up to 
1050nm). The PL intensity is spin dependent, in case of 
Vsi(V2) it is higher when the system is initialized into 
the ms = ±3/2 states. The radiofrequency (RF) pro¬ 
vided by a signal generator (2 dBm power) is guided to 
a thin copper wire and terminated with 50 Vl impedance. 
The laser beam is focused close to the wire using a lOx 
optical objective (N.A. = 0.25) and the PL is collected 
through the sample using a biconvex lens. The x-axis 
of the laboratory coordinate system is set parallel to the 
wire. The external magnetic field B can be applied along 
arbitrary directions, using a 3D coil arrangement. The 
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FIG. 1: (a) The c-axis of the 4H-SiC crystal is oriented at 

an angle of 7° with respect to the surface normal. The sili¬ 
con vacancy Vsi(V2) has spin-3/2 ground state with zero-held 
splitting uq = 70 MHz. (b) Room-temperature ODMR spec¬ 
tra of the Vsi(V2) defect in different magnetic helds applied 
parallel to the surface normal {B\\z). 


magnet is calibrated using a 3D Hall sensor, providing 
angle resolution of 2 °. 

Without external magnetic held, a resonance RF z/q = 
2D/h induces magnetic dipole transitions between the 
spin-split sublevels (±1/2 ^ ±3/2) [Fig. [^a)], result¬ 
ing in a change of the PL intensity (APL). To increase 
the sensitivity, the RF is chopped at 677 Hz and the out¬ 
put PL signal is locked-in. An example of the ODMR 
spectrum (i.e., the ODMR contrast APL/PL against ap¬ 
plied RF frequency) obtained for 5 —> 0 is shown in 
Fig-Bb) (lower curve). The ODMR line is split around 
z/q = 70 MHz due to the geomagnetic held and stray mag¬ 
netic helds. 

We now discuss the evolution of the ODMR spectra in 
external magnetic helds. The corresponding spin Hamil¬ 
tonian is written in the form 

U = ge^B'^S + D[Sl-S{S + l)/^]. ( 1 ) 

Here, ~ 2.0 is the electron g-factor, is the Bohr 
magneton and Sc is the projection of the total spin on the 
c-axis of 4H-SiC [27]. For the sake of simplicity, we do 
not consider the hyperhne interaction. We also neglect 
any deviation from the uniaxial symmetry, described by 
the transverse zero-held splitting parameter E D. Re¬ 
markably, the ms = ±1/2 {ms = ±3/2) states remain 
doubly degenerated even in the presence of electric and 
strain helds, in accord with the Kramers theorem. We 
hence take £■ = 0 as a good approximation Eilisi]. 

When external magnetic held is applied parallel to the 
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FIG. 2: Spin sublevels in the Vsi(V2) ground state, cal¬ 
culated using Hamiltonian Q as a function of the magnetic 
held B applied (a) parallel to the c-axis and (b) at an angle 
0 = 20° with respect to the c-axis. Magnetic held evolution 
of the ODMR spectrum calculated for (c) ^ = 0° and (d) 
0 = 20°. The ODMR contrast is coded by the line thickness. 
The transitions, associated with the m ODMR line, are shown 
by arrows in (a) and (b). 


symmetry axis B\\c {0 = 0°), the spin states are split as 
e±i /2 = -D ± geHBB/2 and €± 3/2 = H--D ± ^geHBB/2 
[Fig. i a)]. One of the dipole-allowed transition with 
Am^ = —1 is (—1/2 ^ “3/2) and the correspond¬ 
ing ODMR line shifts linearly with magnetic held z/i = 
1^0 “ Qel^B^/h\^ as shown in Fig. [^c). Another dipole- 
allowed transition with Am^ = +1 is (±1/2 ^ ±3/2) 
and the corresponding ODMR line shifts linearly towards 
higher frequencies V 2 = r'o ± del^BB/h. For magnetic 
helds lower than Bq = hvo/{ qcIIb) = 2.5 mT the split¬ 
ting between the ODMR lines yields z ^2 = ‘^Qe^BBjh 

[middle curve in Fig.j^b)] and for higher magnetic helds 
B > Bq this splitting is equal z ^2 — M = 4D//z [upper 
curve in Fig.j^b)]. Remarkably, the dipole-allowed tran¬ 
sition (— 1/2 + 1 / 2 ), expected between ui and z /2 at 

frequencies gd^BB/h^ is not observed in the ODMR spec¬ 
tra. The reason is the equal population of the ms = ±1/2 
sublevels due to the optical pumping. 

The behavior becomes much more complex when the 
magnetic held vector B is no more parallel to the defect 
symmetry axis. It is instructive to consider a certain pro¬ 
nounced case in detail, namely when the magnetic held 
is applied at an angle 0 = 20 ° with respect to the c-axis. 
By solving Hamiltonian 0 , we hnd the energy level posi- 























3 



0 20 40 60 80 100 120 140 160 180 200 


Frequency, v (MHz) 

FIG. 3: Evolution of the Vsi(V2) ODMR spectra in magnetic 
fields applied parallel to the surface normal i.e., at an 

angle 0 — 7° with respect to the c-axis: (a) experiment and 
(b) calculation. The ODMR contrast is color-coded (bright 
colors correspond to higher values). 


tions and the eigenfunctions \j) (j = 1, 2, 3,4) depending 
on B [Fig. [^b)]. The mixing of the ms = ±1/2 states 
in the transverse component of the magnetic field B sin 0 
results in the appearance of two additional ODMR lines, 
which are labeled z /3 and z /4 in Fig. [^d). The level anti¬ 
crossing seen in Fig. |^b) at 1.3 mT (between states |2) 
and |3)) and at 2.0 mT (between states |1) and |2)) re¬ 
veals as the appearance of ’turning points’ at frequencies 
j^min _ and = 39 MHz in the calculated 

spectra of Fig. id). Furthermore, according to Fig. [^b) 
the population of states | 1 ) and | 2 ) becomes equal at 
3.0 mT and the ODMR contrast tends to zero, seeing as 
a discontinuity at = 47 MHz in Fig. |^d). 

We now explain how the field-dependent ODMR spec¬ 
tra in Fig.j^d) are calculated. To obtain the probability 
of the RF-induced transitions between states \j) and |/c), 
we apply 


Wjk ^\{j\B^S^\k)\\ (2) 

Here, Bi is the driving RF-field (in our experiments 
We modify the spin-3/2 matrices Si in Eq. Q to 
take into account that the ms = ± 1/2 states are equally 
populated due to the optical pumping and hence do not 
contribute to the ODMR signal. Namely, the matrix el¬ 
ements coupling these states are set to zero. Finally we 
simulate the ODMR spectra assuming Lorentzian line 
shapes with the experimental value for the full width 
at half maximum of 2.7 MHz m- One can see from 
Fig.|^d) that for an arbitrary strength and orientation of 
the magnetic field up to six ODMR lines can be observed. 

Using Eq. (§ with modified spin matrices we can 
perfectly describe our experimental data. Figure ia) 
presents the evolution of the ODMR spectra when the 



FIG. 4: (a) Evolution of the Vsi(V2) ODMR spectra in 

magnetic fields applied at an angle 0' = 60° with respect to 
the surface normal in the yz-plane: (b) experiment and (c) 
calculation. 


external magnetic field is applied parallel to the z-axis of 
the laboratory coordinate system, i.e. at an angle 0 = 7° 
with respect to the c-axis [see Fig. [^a)]. Because the 
deviation from the defect symmetry axis is small, the be¬ 
havior is very similar to the high-symmetry case {B\\c) 
of Fig. [^c). The difference is the manifestation of anti¬ 
crossing = 13 MHz) between the {ms = —3/2)-like 

and {ms = —l/2)-like states at 2.5 mT. Furthermore, in 
addition to the inner resonances and 1 ^ 2 , the outer res¬ 
onances Us and z /4 appear. They shift with double slope 
of the inner resonances. This behavior is closely repro¬ 
duced by our calculations shown in Fig. [^b). 

As another demonstration, we measure the evolution 
of the ODMR spectra when the external magnetic field is 
applied in the ^z-plane at an angle 0 ' = 60° with respect 
to the 2 )-axis [Fig. [^a)]. This corresponds to the angle 
0 ^ 61° between B and the c-axis. It is close to the magic 
angle Om = arccos(l/\/3) ~ 54.7°, when the splitting be¬ 
tween the inner ODMR resonances vanishes. Indeed, we 
experimentally observed U 2 — ui 1 ^ 4 ^ — us [Fig. |^b)]. 
Furthermore, when the magnetic field is significantly in¬ 
clined from the symmetry axis, the outer resonances z /3 
and z /4 become much more pronounced. This compli¬ 
cated behaviour is reproduced in the calculated ODMR 
spectra of Fig. |^c) as well. 

We now discuss how the orientation of the magnetic 
field can be reconstructed. Figure [^a) shows the evo¬ 
lution of the ODMR spectra depending on the magnetic 
field orientation at a fixed strength H = 0.5 mT. The 
orientation angle varies from 0' = 0° (B||z) to 0' = 90° 
(B± 2 )) in the x 2 )-plane. As the defect symmetry axis is 
not exactly parallel to the z-ax.is and lies in the x 2 )-plane 
[the schematic of Fig. [^a)] we shift the calculated re¬ 
sults accordingly {0 = 0' — 7°) for the direct comparison 
[Fig.§b)]. As one can see from Figs, [^a) and (b), the 
splitting between the inner { 1^2 — ui) and outer (z /4 — z/ 3 ) 
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FIG. 5: Evolution of the Vsi(V2) ODMR spectra in a mag¬ 
netic field B = 0.5 mT revolving in the xz-plane: (a) experi¬ 
ment and (b) calculation. The vertical axis in (b) is shifted by 
7° to account for the c-axis tilting with respect to the z-axis. 
(c) Polar angle dependence of the relative splittings between 
inner (z /2 — J^i) and outer (z /4 — 1 ^ 3 ) resonances. The dashed 
line corresponds to the angles when the outer ODMR lines 
are weak. 


resonances changes with magnetic field orientation. Fig¬ 
ure [^c) represents the interconnection between the polar 
angle 6>, i.e., the angle between the magnetic field vector 
and the c-axis, and the relative splitting between reso¬ 
nances tz = ~ ^ 3 )- 

Thus, the algorithm to determine the orientation of 
the magnetic field can be described as follows. For the 
relatively large polar angles 0 > 54.7° the ms = ±1/2 
states are strongly mixed and the outer resonances are 
well resolved. To determine 0 in this case, one measures 
the relative distance between the ODMR lines k and uses 
the calibration curve shown by the solid line in Fig. [^c). 
In the opposite case when 0 < 54.7°, the ms = ±1/2 
states are weakly mixed and the outer ODMR lines have 
lower amplitude, less than 50% compared to the inner 
ODMR lines. The corresponding calibration curve to de¬ 
termine 0 is shown by the dashed line in Fig. gc). The 
angle resolution depends on the distance between the spin 
resonance frequencies at a given magnetic filed and the 
ODMR linewidth. For the data presented in Fig. the 
anglular resolution is better than 5°. It should be pos¬ 
sible to improve it further by using advanced readout 
protocols [2]. Remarkably, at 6> = 0° and 0 = 90° only 
two ODMR lines should be observed and therefore they 
are formally indistinguishable. However, in the presence 
of weak stray fields they behave differently. For 0^0° 
the outer resonances with double splitting become visible 


[as in Fig. |^a)] and, for 0 90°, the outer resonances 

are of the same amplitude as the inner resonances and 
only slightly split from them. 

In summary, uniaxial spin-3/2 centers in hexagonal lat¬ 
tice can be used to measure magnetic field strength and 
to reconstruct magnetic field orientation with respect to 
the symmetry axis of the crystal (i.e., the polar angle). 
The method is based on the ODMR technique and ap¬ 
plicable for ensembles of spin centers as well as for single 
centers. As a probe, we used the Vsi spin center in 4H- 
SiC and performed demonstrative experiments in weak 
magnetic fields up to 5mT. The experiments are per¬ 
fectly described by our model for any magnetic field ori¬ 
entation. It should be mentioned that this method does 
not provide information on the azimuthal angle, but this 
can potentially be overcome by using differently oriented 
centers, for instance along the c-axis and in the basal 
plane of the crystal. Interestingly, Vsi can be incorpo¬ 
rated into SiC nanocrystals [33] and their density can be 
controlled over eight orders of magnitude down to sin¬ 
gle defect level [32]. The selective optical addressability 
m and coherent control [12 of single Vsi centers m 
have also been reported. Combination of these capabil¬ 
ities with our findings suggests promising perspectives 
for vector magnetometry and local imaging down to the 
nanoscale. 
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